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Abstract 

Computation of the renormalized stress-energy tensor is the most serious obstacle in studying 
the dynamical, self-consistent, semiclassical evaporation of a black hole in 4D. The difficulty arises 
from the delicate regularization procedure for the stress-energy tensor, combined with the fact 
that in practice the modes of the held need be computed numerically. We have developed a new 
method for numerical implementation of the point-splitting regularization in 4D, applicable to the 
renormalized stress-energy tensor as well as to {f 2 ) re , namely the renormalized (f 2 )- So far we 
have formulated two variants of this method: t-splitting (aimed for stationary backgrounds) and 
angular splitting (for spherically-symmetric backgrounds). In this paper we introduce our basic 
approach, and then focus on the t-splitting variant, which is the simplest of the two (deferring 
the angular-splitting variant to a forthcoming paper). We then use this variant, as a first stage, 
to calculate (f 2 ) ren hi Schwarzschild spacetime, for a massless scalar held in the Boulware state. 
We compare our results to previous ones, obtained by a different method, and find full agreement. 
We discuss how this approach can be applied (using the angular-splitting variant) to analyze the 
dynamical self-consistent evaporation of black holes. 
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I. INTRODUCTION 


After the discovery of Hawking radiation [Tj in 1975, it was widely anticipated that the 
semi-classical approach to gravity based on quantum field theory (QFT) in curved space- 
time will open the opportunity to explore various interesting physical problems, in which a 
quantum held interacts with the curved spacetime where it resides. Among these problems, 
of special interest is the self-consistent dynamical evaporation of a black hole (BH). Another 
outstanding problem in this class is the evolution of quantum fluctuations in the very early 
universe, and the resulting cosmological structure formation. 

Indeed, in the latter problem of cosmological quantum-held perturbations a remarkable 
progress has been achieved over the last few decades. M However, the problem of analyzing 
the semiclassical evolution of an evaporating BH still remains a serious challenge. 

To understand the difficulties in analyzing this interesting problem of self-consistent BH 
evaporation, let us briehy review the basic structure of semiclassical gravity. The metric 
g a p{x) is treated as a classical held, and it is assumed to satisfy the semiclassical Einstein 
equation 

G a/3 = 8vr (T a p) ren , (1.1) 

where G a/ 3 is the Einstein tensor associated with g a / 3 (x). Q The source term (T a p) ren is the 
regularized expectation value of the stress-energy tensor associated with a quantum held 
(pfx). For the sake of simplicity we shall take (p here to be a scalar held. It is supposed to 
satisfy the held equation 

(□- m 1 2 -^)0 = 0 , ( 1 . 2 ) 


where m and £ respectively denote the mass and coupling constant of the scalar held. Note 
that the theory is semi-classical, as the held (p is quantized but the metric g a p is classical, 
nevertheless the two are coupled. 

The major difficulty in solving (or even analyzing) the held equation (1.1) has to do with 
the regularization of the divergent quantity (T a p). [^] As opposed to QFT in hat spacetime in 
which one can use the normal-ordering procedure, in curved spacetime the outcome of this 
procedure depends on the choice of time slicing which is completely arbitrary. There exists a 
regularization method named point-splitting (PS), also known as covariant point separation, 
which gives a general prescription how to regularize quantities which are quadratic in the 
held and its derivatives such as (T a p). Alas, implementing this prescription in situations 
where the solution of the held equation (1.2) is known only numerically turns out as a 
surprisingly difficult problem. 

In this paper we shall focus on the regularization of ((p 2 ) instead of (T a p). The quantity 
(0 2 ) is also divergent, but not as strong as (T a p). In addition, the scalar character of (<p 2 ) 
(as opposed to the tensorial (T a p)) makes it easier to regularize. These properties make (cp 2 ) 
a convenient tool to examine and explain new ideas concerning regularization. In order to 
calculate (cp 2 (x)) ren (i.e. the regularized (cp 2 (x))) using the PS method we split the point x 
and write it as a product of (p at two different points, namely (cp (x) (p (x')). This is known as 
the two-point function (TPF). We then subtract from the TPF a known counter-term and 
take the limit x' —> x. This limit, however, is what makes the numerical implementation so 
hard. 

In 1984 Candelas and Howard [7j developed a method to numerically implement PS if 
a high-order WKB approximation for the held modes is known. Using this method, they 


1 Throughout this paper we use relativistic units c = G = 1 and the (—b ++) signature. 

2 There are also other difficulties, e.g. the runaway problem, as noted by Wald [4]. 
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calculated {4> 2 ) ren in Schwarzschild spacetime, and subsequently this method was used by 
Howard [8] to calculate {T a p) in Schwarzschild. Later, this method was extended to 
a general static spherically symmetric background, first for (4> 2 ) ren by Anderson |JJ, and 
subsequently to (T a p) ren by Anderson, Hiscock and Samuel [10]. 

The limitation of the method proposed by Candelas and Howard is that it requires a 
high-order WKB approximation for the held modes: at least second-order for (4> 2 ) ren and 
forth order for (T a p) . In the ordinary (Lorentzian-signature) Schwarzschild metric this 
task of high-order WKB expansion is very difficult, especially because of the presence of a 
turning point: For typical modes of large c o and l , there is a turning point on the r axis, 
at a value r turn (u,l) where the effective potential Vfr) equals u 2 . [^] The mode’s radial 
function is essentially oscillatory at r > r turn and exponential at r < r tU m- Both these basic 
WKB approximations — the oscillatory approximation at r > r turn and the exponential 
approximation at r < r turn — break down and actually diverge at r —y r turn . To correctly 
match the two approximations, one has to use another, intermediate approximation valid 
in the neighborhood of r = r turn . This turning-point approximation is based on the Airy 
function. Whereas the leading-order matching is manageable, it becomes exceedingly hard 
to go to higher-order WKB, because each succeeding order will now require its own turning- 
point matching. Furthermore, to the best of our understanding, the series of powers involved 
in the Airy-based turning-point expansion proceeds in powers of u ; -1 / 3 (rather than a; -1 ). 
Correspondingly, we may expect that in the presence of a turning point, to implement the 
WKB-based expansion to order u ; -4 (required for calculating (T a p) ), one would have to 
carry matched asymptotic expansion up to twelfth order in cu -1 / 3 (or six such orders for 
(4> 2 ) ren ) — a formidably difficult task. 

To overcome these difficulties, Candelas and Howard [7] and several others [EHIQJ used an 
elegant trick: They used Wick rotation to analytically extend the background metric to the 
Euclidean sector. Any static spacetime is guaranteed to have such a real-metric Euclidean 
sector. In the latter, the radial equation does not admit a turning point. This way, it 
was possible to carry the WKB analysis to the desired order and to implement the above 
regularization scheme — for ( 4> 2 ) ren as well as for {T a p) ren . 

Our ultimate goal, however, is to develop a regularization scheme for (f 2 ) and (T a p), 
applicable to time-dependent backgrounds as well. Such a time-dependent metric (even if 
spherically symmetric) does not generically admit a Euclidean sector. We are thus led to 
carry the analysis directly in the Lorentzian sector, which in turn implies the presence of 
a turning point in the radial equation, hampering any attempt to carry high-order WKB 
expansion. We shall therefore refrain from establishing our regularization scheme on the 
WKB analysis. 

There is another obvious reason for avoiding WKB analysis: Consider a time-dependent 
spherically-symmetric background. The held equation for a given mode may still be ex¬ 
pressed as a one-dimensional (namely 1 + 1) wave equation with an effective potential, but 
now the potential will be time-dependent. In such a situation, even the leading-order WKB 
(and even if we forget for the moment about the turning point) becomes a non-trivial task, 
let alone higher-order WKB analysis. 

For these reasons, we shall not base our regularization scheme on high-order WKB ex¬ 
pansion. Instead, in our method we extract the required information concerning the liigh- 

3 In Schwarzschild spacetime (and any other eternal BH spacetimes), for given ui and l there are usually 
two such roots. Here we shall explicitly refer to the larger one, but the same complications arise also at 
the smaller root. 
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frequency field’s modes directly from the well-known counter-term ( |2.6 ) for (0 2 ) ren (and, for 
computing ( T a g) r , from the counter-term developed by Christensen 6]). 

We point out that this approach, namely extraction of the high-frequency asymptotic 
behavior of the modes from a known local counter-term, was actually initiated by Candelas 
HU - already before he and Howard resorted to the Euclidean sector [7]. However, Cande¬ 
las’ analysis was restricted to the Schwarzschild case (which in particular means restriction 
to staticity, spherical symmetry, and to vacuum^]). We should also comment that even in 
that case the analysis in Ref. HU was not completed, because the required integral of the 
regularized mode contribution over uj was not carried out. When we attempted to imple¬ 
ment this integral over uj, we found that it actually fails to converge in the usual sense, 
due to growing oscillations (see below), which led us to introduce the notion of generalized 
integral. Our approach is in this sense a completion of Candelas’ method, as well as its gen¬ 
eralization beyond the vacuum case (and with the scope of further extending it to dynamical 
backgrounds). 

Our method requires the held modes to admit a trivial decomposition in at least one of 
the coordinates (e.g. through e~' lujt or spherical harmonics). This usually corresponds to 
having a Killing held in spacetime. [)] The splitting is then done in that trivial coordinate — 
which enables us to treat the coincidence limit analytically. So far we have developed two 
different variants of our method: (i) the t-splitting variant, which requires a time-translation 
symmetry; (ii) the angular-splitting variant, which requires spherical symmetry. We are also 
exploring a third variant, azimuthal splitting (which would only require axial symmetry), 
but this one is still in progress. In all these variants we assume for simplicity that the 
background is asymptotically-hat, although this requirement can probably be relaxed. 

Since our ultimate goal is to analyze ( T a p) ren on the time-dependent background of an 
evaporating BH, the t-splitting method is insufficient, and we shall actually need the angular- 
splitting variant. It turns out, however, that the t-splitting variant is in some sense concep¬ 
tually simpler and easier to present at first stage, because certain additional complications 
arise in the angular-splitting method. Although we know how to address these complica¬ 
tions, they make the method’s logical structure a bit more obscure and harder to explain. 
For this reason, for the sake of introducing our basic regularization strategy we choose to 
present here the t-splitting variant, which is logically simpler. And for exactly the same 
reason, in this paper we shall display the regularization of (0 2 ) rather than {T a g). Then in 
the next paper we plan to present the angular-splitting variant, which is going to be our 
main tool for analyzing dynamical BH evaporation. 

The TPF diverges for any pair of points connected by a null geodesic, even if they are 
far from each other [12j. As it turns out, this long-distance divergence of the TPF leads 
to undamped oscillations in the mode contributions at large uj. To address this issue we 
use the concept of generalized integral , in which these oscillations are properly damped upon 
integration over uj, which fully cures the oscillations problem. The origin of this complication 


(the presence of connecting null geodesics) is discussed in Sec. IIB and also in Appendix 
|B| And the resolution of the oscillations problem by means of generalized integral (and 
particularly the so called “self-cancellation integral”) is described in Sec. IIB and further in 
Appendix [A} 

In the description of the t-splitting method in Sec. 


Ill we assume a spherically-symmetric 


static background for the sake of simplicity. However, as was discussed above, the t-splitting 

4 In particular, the logarithmic counter-term is not encountered in the vacuum case. 

5 Having at least one trivial coordinate is a necessary condition for the applicability of our method, but we 
do not claim that it is also a sufficient condition. 
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method does not require spherical symmetry, and in principle it may be applied to a generic 
(asymptotically-flat) stationary spacetime. We outline this generalization of the method to 
stationary backgrounds in Sec. Ill B We point out, however, that some completion is still 
required in the case of a stationary eternal BH (see therein). 

Next we apply our method explicitly to the Schwarzschild case, computing (4> 2 ) ren in 
Boulware state. We compare our results to those obtained previously by Anderson (using 
the Euclidean sector), and find full agreement. 

This article is divided as follows: In Sec. [IT] we present the basic PS method, and then 
briefly outline the procedure developed by Candelas and Howard [7]- Note that Sec. IIB 


discusses certain subtleties of the TPF in some detail, and may be skipped in Erst reading. 
In Sec. |III| we present our /-split t ing method. We first describe it for spherically-symmetric 
static background, and then outline its generalization to a generic stationary background. 
In Sec. |TV we harness this method for the calculation of ( 4> 2 ) ren hi the Schwarzschild metric. 


Finally, in Sec. [V]we discuss the implications of our new method and try to pave the path 
towards our ultimate goal of investigating self-consistent BH evaporation. 


II. BASIC POINT-SPLITTING METHOD AND ITS NUMERICAL IMPLEMEN¬ 
TATION 


We start by sketching the basic PS regularization method. This method is aimed to 
regularize the expectation value of various quantities which are quadratic in the held operator 
(and its derivatives). Among these quantities, the most important one is probably the 
energy-momentum tensor ( T a p ). However, in this first paper we shall consider (0 2 ) as a 
simpler example (although we shall occasionally remark on the analogous calculation of 

For simplicity we consider here a quantum scalar held 4>(x) living in a static, spherically 
symmetric, asymptotically hat spacetime with metric 

ds 2 = g tt (r)dt 2 + g rr {r)dr 2 + r 2 dVt 2 , (2.1) 

where dVt 2 = dd 2 + sin 2 6 dtp 2 . The held operator may then conveniently be expressed as 


(p(x) 



OO l 

duj ££( fculm (^) &L dim “t - fujlm 

1=0 rn=—l 



( 2 . 2 ) 


Here, a) ulm and a^im are the creation and annihilation operators of the field’s culm mode, 
fujlm {%) is a complete, orthonormal, family of modes [^] taking the form 


fulm (®) = e Y im ( e , (p) $ ul (r), 


(2.3) 


and Yi m (0,(p) are the usual spherical harmonics. The radial functions 'i/y, i (r) are obtained 
by solving the held equation for 4>(x) with the decomposition (2.3) and with appropriate 
boundary conditions. 

The operation of summation over m, l and integration over uj repeats many times in the 
analysis below. We shall generally refer to this operation as the “mode sum” (despite the 


6 By “orthonormal” we mean that the inner product of two mode functions and is Su>S m m'S(ui — 
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slight abuse of terminology). We point out, however, that the decomposition (2.2) applies 
as-is in the case of an asymptotically-flat background spacetime with simple asymptotic 
structure (like e.g. Minkowski or a star). But if the background spacetime is an eternal 
BH with a past horizon, then for each a dm combination there are actually two orthonormal 
modes, namely the “in” and “up” modes. The “in” modes are those described above (namely 
monochromatic waves propagating from past null infinity), and the “up” modes describe 
monochromatic waves that emerge from the past horizon. In this case of eternal BH, the 
mode sum should also include a summation over the contributions of these two independent 
modes for each ulm (as explicitly described in Sec. IV for the Schwarzschild case). 

The quantity (</> 2 ) obviously depends on the quantum state. Naturally, one would like to 
evaluate it in the vacuum state. The latter is defined to be the quantum state annihilated 
by each of the above a u i m operators. Q 

Trying to naively calculate (</> 2 (x)) in the vacuum state yields the divergent expression 


(0 2 (a)) r 


= h 


, OO l 

\ Ylm (^’^l 2 l^/i ( r )T 

1=0 rn=—l 


(2.4) 


Although the sum over m, l does converge for a given u>, the integral over u diverges. In fact, 
one can easily check that already in Minkowski spacetime the integrand is oc a; (and oc cu 3 
for (T a p)), and the same divergence occurs in the Schwarzschild case as well. If one tries to 
integrate over u ; before the summation, one finds that the integral over u> again diverges (for 
a given l, m). 

One therefore needs to somehow regularize the expression for {(ft 2 ). In flat spacetime this 
regularization may be easily achieved by normal ordering, but in curved spacetime this yields 
a slicing-dependent, non-unique result. This is where point splitting comes into play. DeWitt 
[5] proposed that (cj) 2 (x)) (and, more generally, quantities quadratic in the held operator and 
its derivatives) can be treated by taking the product of the held operators in two separate 
points x,x' and then considering the coincidence limit x —> x' . More specifically, he showed 
that the regularized expectation value of </> 2 can be dehned as 

( X ))ren = [(0 0) 0 0')) “ G DS (x, x')} . (2.5) 

Here Gds (x,x') is the DeWitt-Schwinger counter-term, , namely a local term which fully 
captures the singular piece of the TPF. For a scalar held with mass m and coupling constant 
£ it takes the form 


= sk 


+ 


m~ + (£ — 1/6) R 
87 r 2 


7 H— In 
[ 2 


/ U 2 (J 


m 


;R 


<j’ a ,h 


167T 2 ' 967r 2 0/3 cr 

Here R and R a p are respectively the Ricci scalar and tensor, 7 denotes the Euler constant, 
and cr (x, x') is the biscalar associated with the short geodesic connecting x and x'. The value 
of a is half the geodesic distance squared (see Ref. m)- More specifically, for a timelike 
separation a = —r 2 /2, where r denotes the proper time between x and x'. The parameter /j 
is unknown and it corresponds to the well-known ambiguity in the regularization procedure 

1.0 


7 If a past horizon exists, then one need to further specify this vacuum state, e.g. by prescribing the outcome 
of the action of the “up” annihilation operators on that state. In the analysis of the Schwarzschild case in 


Sec. IV we shall consider the Boulware vacuum state. 

8 In the case of a massive scalar field some authors took /i to be the field’s mass m. Gann] 
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A. Previous numerical implementations of point splitting 


It is not easy to directly implement the point-splitting method, especially because in 
the cases of interest (e.g. black hole backgrounds) the radial functions are only known 
from numerics. Furthermore, in Eq. (2.5) the TPF ((ft (x) (ft (x')), which is to be computed 
numerically, diverges as r~ 2 (and as r^ 4 for (T a p)) as x' approaches x , where r denotes the 
geodesic distance between x and x'. 

This problem was addressed by Candelas, Howard and later Anderson and collaborators 
(see Refs. HHH3]) a long time ago. They developed a calculation scheme that allows regu¬ 
larization of ((ft 2 ) or (T al s) numerically, provided that an analytic approximation for the held 
is known, up to a sufficiently high order [^J 

In that scheme, one first analytically constructs the approximate singular piece of the 
held, which we denote <ft S ing(x). It is composed of the contribution of the modes of large 
uj and l, up to second order in 1/co and l/l (fourth order for (T a p)), which is usually 
computed using WKB analysis. From this quantity one then constructs the approximate 
TPF ((ft singsing(x')), and recasts Eq. (2.5) in the form 


{(ft 2 ( X ))ren = J)™ 1(0 ( X ) 0 0')) “ ((ftsing(x)(ft sin g(x'))\ 

+ hill [{(ftsing(x)<ftsing(x')) ~ G DS (x,x')\ . (2.7) 

rp! -V rp 


Since (<; ft sing (x) (ft sing (x ')) contains the entire singular piece of the TPF, both limits at the 
R.H.S. are well defined. We denote the second limit by (eft 2 (x)) Anolyllc . This quantity is 
well-defined and regular, and it is computed analytically by summing/integrating over the 
WKB expressions for the large oj, l modes which comprise (ft S i ng . The hrst term in the R.H.S. 
is now expressed as a mode sum, and owing to its regularity the limit may trivially be taken 
by replacing x' by x. Putting it all together one obtains 


/»00 OO ^ 

(<l )2 W)ren = h dUJ J2 £ l^m (0, ifi) | 2 (| $ ul ( t ) f ~ (r) f) + {(ft 2 (x) > . 

Jo 1=0 m=—l 

( 2 . 8 ) 

The mode sum should now converge (even though it is taken in coincidence), owing to the 
subtraction of the singular piece. 

Using this method with splitting in the t direction, (<ft 2 ) ren (as well as ( T a p ) ren ) was 
calculated for the Schwarzschild case HI El , and later also for a generic static spherically- 
symmetric spacetime mm- 

Besides the need to numerically compute the various mode functions (r), this method 
also includes a challenging analytical component: It requires a high-order WKB analysis. As 
was mentioned in the Introduction, the presence of a turning point makes this an extremely 
difficult task. To overcome these difficulties, Candelas and Howard [7j and others PEoj 
actually carried the analysis in the Euclidean sector. 


9 In general it would require a second order WKB approximation in order to compute < cf i 2 > re n, and a 
fourth order approximation to compute <T a p > re n- 
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B. Remarks about the TPF and its mode sum F^I 

1. Caveats concerning the regularity of the TPF itself 

The point-splitting method is based on the presumption that the two-point function is 
well-behaved as long as the two points are separated. A few caveats are associated to this 
issue: First, even in flat spacetime, the two-point function diverges when the separation 
is in a null direction. Let us therefore assume, for the sake of simplicity, that the points 
are separated in a timelike (or possibly spacelike) direction. The second caveat is that in 
curved spacetime, assuming that the two points are indeed separated by a timelike geodesic, 
if the (proper-time) distance between x and x' is sufficiently large, there may also be a null 
geodesic connecting these two points. For example, consider the Schwarzschild spacetime 
and an approximately-static timelike geodesic I located in the asymptotic region very far 
from the BH. Let i be a point on I . There is a null geodesic which emanates from x and 
moves towards the BH (but with an appropriate miss), makes a turn around the BH, then 
returns to T and hits it at a point X\. It turns out that (0 ( x ) 0 (x')) develops a singularity 
at x' = xi D2P- In fact, along the timelike geodesic T there is an infinite discrete set of 
points x n which are connected to x by null geodesics that make n turns around the BH 
before returning to T, and (0 (x) 0 ( x ')) diverges at all points x' = x n . 

Note, however, that the problematic points x n are all located far away from x, outside the 
normal neighborhood; and for the point-splitting procedure only points x 1 in the immediate 
neighborhood of x are relevant. We shall thus restrict our attention now to points x’ in the 
close neighborhood of x, with a timelike (or alternatively spacelike) separation between x 
and x'. Then the TPF should be well-behaved. 


2. Caveats concerning the convergence of the mode sum 


Naively one might expect that since the TPF is well behaved for a short timelike or 
spacelike separation, the mode sum associated to it should converge. It turns out, however, 
that the situation is more subtle: The involved sum/integral usually fails to converge (in the 
literal sense). Typically this failure to converge is associated with undamped oscillations. 

This phenomenon, the failure of strict convergence, occurs already in flat spacetime. 
The (non-) convergence situation may depend on the direction of splitting, on the order of 
operations (summation over l , m and integration over cu), and on the specific quantity being 
calculated (whether it is the TPF, or a component of ( T a p )). As an example, consider the 
calculation of (0 (x) 0 (x')) in Minkowski spacetime, using standard mode decomposition in 
spherical harmonics and temporal modes e~ lujt . The mode functions are given by 


0 ui (r) 


Wr J,+1/2 K) ’ 


where J is the Bessel function of the first kind. Consider now a point x located at some 
r > 0, and a point x' displaced in the t direction by an amount t' — t — At. The sum over 
Z, m (for a given u) then converges, and one finds that 


i 


OO 

£ 

1=0 m=—l 


\Yim(0,T)\ 2 \^ui(r)\ 2 =—u 
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This subsection is somewhat remote from the main line of this paper, and can be skipped at first reading. 




(2.9) 


The mode-sum expression for the TPF then becomes 

t coo 

{(/) (x) (j) 0 ')) = — J coe luAt dcj. 

The problem is that this integral does not converge in the usual sense (even conditionally), 
due to the growing oscillations at large co. 

In this example of pure f-splitting, in which we first sum over m, l, it was the integral over 
co that failed to converge. In another application of point splitting, in which the splitting is 
in both t and 6, and in which one first integrates over co and sums over m, l afterward, one 
finds that this time the sum over l fails to converge, again due to growing oscillations. [^] 
This situation, of non-converging oscillatory integrals (or sums over /), is fairly common 
in various QFT calculations. The common practice (which may be justified by several 
arguments) is to contend that the large-a; oscillations should be damped in some appropriate 
manner. 


3. Generalized integral 


This situation motivates us to introduce the notion of generalized integral, which properly 
incorporates oscillation damping. We should emphasize that all the co integrals in this paper 
are in fact generalized integrals. 

We shall consider here two specific procedures of oscillation damping, which yield two 
(mutually consistent) definitions of generalized integral. 

Abel-summation integral: This is a commonly-used method for giving a meaning for such 
oscillatory integrals. The Abel-summed integral is defined as 


r- 00 (A) 


'0 


h(co)dco = lim / e eoJ h(co)dw 

e_5>0 + Jo 


( 2 . 10 ) 


- provided of course that the integral at the R.H.S. is well defined for e > 0, and the limit 
e —» O 4 . exists. For example, with Abel summation, the integral in Eq. (2.9) reads — 1/At 2 


- yielding the standard expression —h/ ( 47 t 2 A t 2 ) for the flat-space TPF [cf. Eq. (2.6)|. 

A crucial property of the Abel-summation integral is that it is consistent with the standard 
integral. Namely, whenever the function h(co) is integrable in the strict sense, its Abel- 
summation integral coincides with the standard integral of this function. 

Self-cancellation integral: For the calculation of (4> 2 ) ren in a black-hole spacetime we 
shall have to carry (generalized) integral of oscillatory functions h{co) that we determine 
numerically. In such a case, this procedure of Abel summation — which must now be im¬ 
plemented numerically — is fairly inconvenient. Furthermore, in practice we only determine 
h(oo) in a restricted range 0 < co < co max , which makes the Abel-summation integral even 
harder to implement. [^] We therefore find it much more convenient to use another concept 
of generalized integral, which we name self-cancellation (of the oscillations). This type of 


11 In addition, the integral over co converges for the TPF but only conditionally, and fails to converge even 
conditionally for T tt , due to growing oscillations. 

12 To this end one would have to generalize the definition of the Abel-summation integral, so as to combine 
the two (non-commuting) limits e —> 0 and co max -A 00 in an appropriate manner. And it turns out that 
the convergence of this generalized Abel integral with increasing cj max is rather slow. 
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generalized integral is applicable whenever the oscillations have well-defined frequencies— 
which is indeed the situation in our problem (see Appendix [B]) . To formulate this concept, 
we first define the integral function 

H(u) = / h(x)dx. (2.11) 

Jo 

The standard integral may then be expressed as 

poo 

/ h(oo)duo = lim H(oo) 

Jo ^°° 

(whenever this limit exists). Instead, our self-cancellation generalized integral is defined as 


r»00(sc) 


h(oo)doo = lim 

CJ—>■ OO 


H{oo) + H(oo + A/2) 


( 2 . 12 ) 


Here, A denotes the “ wave -length 11 of the oscillation in h{uo) (which is also inherited by H). 
For example, in Eq. (2.9) the oscillatory factor is e lU]At , hence the period of oscillation is 
A = 2n/At. The idea is simple: The non-oscillating piece of H(uo) is unaffected by this 
averaging, but the oscillatory piece will be very effectively annihilated by such averaging 
with half-wavelength shift. 

As a simplest example, consider the case h(oo) = sin (ooL). Then H{oo) = [1 — cos((nL)] /L, 
and obviously A = 27 t/L. Clearly H(oo) fails to have a limit uj —> oo. Yet the self-cancellation 
integral is perfectly well-defined: The term in squared brackets in Eq. (2.12) is simply 1/L, 
entirely independent of oo. This example demonstrates the potential of the self-cancellation 
method to yield extremely fast convergence in oo. This last property is important, especially 
because in an actual calculation we have to determine h(oo) numerically, and we do so in a 
restricted range of oo. 

For later convenience we re-formulate this notion of self-cancellation integral as follows: 


^oo(sc) 


h(uo)duo = lim T\ [H(oo )], 


(2.13) 


where T\ is the “self-cancellation operation” defined by 


Tx [/M] = 


f( u, ) + f(u> + A/2) 


In the actual calculation of (0 2 ) we shall have to repeat the self-cancellation operation 
several times (for several different oscillation frequencies). In Appendix [A] we shall introduce 
the “multiple self-cancellation operation” T*, formed by combining several T\ operations. 

It is easy to show that the notion of self-cancellation integral is fully consistent with the 
standard integral — whenever the latter is well defined. Furthermore, the self-cancellation 
integral is also fully consistent with the Abel-summation integral, in the following sense: If 
the self-cancellation integral converges, then the Abel integral converges too, and the two 
generalized integrals yield the same result. Note, however, that the Abel summation method 
is more general than self-cancellation. Namely, there are functions h{oo) for which the Abel- 
sunnned integral is well defined but the self-cancellation integral is not (that is, the limit 
in Eq. (2.12) is non-existent). Nevertheless, for the functions h(oo) involved in the analysis 


below, the self-cancellation integral is well-defined and extremely powerful. 
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III. OUR NEW METHOD: THE ^-SPLITTING VARIANT 


As was already noted in the Introduction, since our ultimate goal is to address dynamical 
background metrics as well, we shall carry the analysis directly in the Lorentzian sector. 
Hence, due to the inevitable presence of a turning point (and also due to the PDE nature of 
the time-dependent mode equation), we shall avoid using WKB expansion in our method. 
Instead we extract the required information about the high-frequency modes of the held 
directly from the counter-term (2.6). 

We shall present here the ^splitting variant, which requires a time-translation Killing 
held. Although this variant should be applicable to a rather generic stationary asymptotically- 
hat background, we shall restrict our attention hrst to the more specihc case of static 
spherically-symmetric background, for the sake of simplicity. Then in Sec. IIIB we outline 
the generalization of the ^splitting formulation beyond spherical symmetry (and beyond 
staticity), to a more generic stationary asymptotically-hat background. 


A. Spherically-symmetric static background 

We split the points in the t direction, namely 

x = (t, r, 6, ip), x' — (t + £, r, 6, ip). (3.1) 

The TPF then takes the form 


(0 


h 


du e 1 ' 


E 

1=0 rn=—l 


K 


I Y tm \ 2 1 tpui (r) 


(3.2) 


As was already mentioned above, the sum over l , m converges, and we denote it by F (c u, r ). 
In fact, the sum of |T)m| 2 over m yields (21 + l)/4(/r, and therefore 

°° o/ | ~i 

= (3.3) 

1=0 


This function is to be computed numerically. The TPF now reduces to 

POO 

(4> (x) (j) (x')) — h / F (u,r) e lU£ du. 


(3.4) 


(Note that in the coincidence limit e —> 0 this integral would diverge. However, the oscilla¬ 
tory factor e lU£ regularizes it.) Equation (2.5) now reads 


(*)}ren = 


£— 


h F(u, r ) e iuJE du - G DS (e) 


(3.5) 


where by Gds (^) we refer to Gds (%, %') with x’ given by Eq. (3.1). Note that for a given x , 
a is uniquely determined by e. By a fairly straightforward Taylor expansion of the geodesic 
equation (and er) in e, one hnds that Gds in Eq. (2.6) takes the general form 


\g ds (x, x') = a (r) e 2 + c (r) 
n 


ITT 


In (e/z) +7 - y 


+ d ( r) + 0(e), 


(3.6) 
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where a(r), c(r), d(r) are certain (real) functions that depend on the background metric and 
the parameters of the held. 13 The explicit form of these functions is not important in the 


present discussion (although it is certainly needed for the actual calculation of (0 2 ) ren )- 
To proceed, we now decompose the e-dependent terms in Gds using the Laplace trans¬ 
form. We have the following identities: 


In (e/i) = — 


e" 2 = 


cu + fl 


ue 


'du , 


'0 


s du+(j- 7 


+ O (e In e). 


Inserting these identities into Eqs. (3.5|3.6 ) we obtain 


k 2 


(x)) = Mini 

v ' ' ren e—>0 


reg 


(oj, r) e lul£ duj — hd{r) 


where 


and 


Fret (w, r) = F (u, r) - F Hng (u,r) , 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


F, im (w, r) = -a (r) u - c (r) 


Consider now the R.H.S. of Eq. (3.9) 


—— ■ (3.H) 

UJ + fJL 

If the integral of F reg (c o, r ) converges (in either the 
strict sense or the generalized sense), we can interchange the limit and integration, and get rid 
of the e —» 0 limit altogether. We should expect the convergence of the integral of F reg (oj, r), 
because the singular piece F sing (a;,r) has already been removed from F (c u,r). In fact, we 
End that in the Schwarzschild case the integral of F reg (oj, r) indeed converges, although 
only in the generalized sense due to oscillations (see next section). We were unable to 
prove the (even generalized) convergence of the integral of F reg (oj,r), but nevertheless since 
this convergence is naturally expected, and since the Schwarzschild example confirms this 
expectation, we shall hereafter assume that this integral indeed converges in the generalized 
sense. E3 

We therefore write our final result as 


( X ))ren = H 


Freg (jd j f t 


where, recall, the integral over u is a generalized one (as defined in Sec 
plementation of this generalized integral is demonstrated in Sec 


doj — hd (r) (3-12) 

2.2). The im- 
lYl for the Schwarzschild 


case. 


B. The general stationary case 

As was already pointed out above, our t-splitting method does not require the background 
metric to be spherically symmetric or even static: It should be applicable to a generic 
stationary asymptotically-flat background. Here we outline this extension. 

13 The term — *7t/2 appears in the brackets because a is negative. 

14 Note that there is no much risk in making such an assumption, because if for a certain background 
metric this assumption turns out to be false, then the attempt to integrate F reg will demonstrate this 
non-convergence right away. 
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On account of asymptotic flatness, we now choose our coordinates (t, r, 9, p) such that 
the large—r asymptotic metric still takes its standard Minkowski form — dt 2 + dr 2 + r 2 dVt 2 . 
p| Due to lack of spherical symmetry, the expression (2.3) for the held modes is now replace 

by 

L dm ( x ) = e~ lut ^ u i m (r, e, p ), 

where ipum (d d, is a set of solutions to the (^-dependent) spatial part of the held equation, 
which is now a PDE (in r, 9, p) rather than ODE. These solutions are required to be regular 
everywhere, and to satisfy the large-r boundary condition 


4>uim (r, 9, <p) = (rV&ruj) L e ( 9 , <p) + (r, 9, <p ), 


r —y oo) 


where denotes the rehected held which is e lujr */r times some function of 9 and p. 

16 1 Once the modes f w i m {x) were defined, the expression (2.2) for the held operator is 


unchanged. 

The TPF now takes the form 


(4> (x) (j) (x')} = h / duj e we Y j y, |-0o dm(r,9,(p) 

- 20 1=0 m=—l 


(3.13) 


It is important to recall that even though the metric is not spherically symmetric, the sums 
over m,! should still converge (for a given oj): Due to asymptotic flatness, at large r there 
will be a centrifugal barrier « i(? + l)/r 2 in the effective potential, just like in Minkowski, 
preventing the penetration of modes with too large l. Thus we can again define 


F (u,r,9, ip) = y y | $ u im(r,d,<p) 

1=0 rn=—l 


(3.14) 


and the TPF still takes the form ( 3.4[ ) [although now with F (cu,r,9,ip)\. The calculation 
now proceeds just as in the spherically symmetric case — except that all the quantities in 


Sec. (Ill A) that were dependent on r only, now depend on 9 and <p as well. The final result 


IS 

poo 



{(f> 2 ( x )) ren = ft F reg (w , D 9 , (p) duj - h d (r, 9 , p) 

Jo 

(3.15) 

with 


Freg (^, T, 9,p) = F ( U , 7", 9, p) ~ F sing (u, T, 9, p) , 

(3.16) 

and 


Fsing (w, r, 9, p) - a{r,9,p)uj c(r,9,p ) . 

uj p 

(3.17) 


It should be pointed out that although the choice of the coordinates r, 9, ip is certainly 
non-unique (due to lack of spherical symmetry), the resultant mode decomposition is still 
unique. This is because of asymptotic flatness (and the standard weak-field metric that 
we require our t,r,9,ip coordinates to satisfy at large r), and because the mode functions 
ipuim (d 9, p>) are defined through their asymptotic form at r —>■ oo. The unique mode 


15 Furthermore we choose our coordinates such that gu takes the standard weak-field form — 1 + 2 M/r + 
0(l/r) 2 , where here M denotes the system’s asymptotic mass. 
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Here r* is to be regarded as the standard function of r given in Eq. (4.41. 
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decomposition in turn leads to a unique Fock space associated to it, and to a well-defined 
vacuum state. 

Note, however, that in the above formulation we implicitly assumed a stationary back¬ 
ground metric with no past horizon (e.g. a spinning star), hence it was sufficient to construct 
the “in” modes. In the case of eternal non-spherical BH the situation becomes more subtle, 
because now we also need to specify the “up” modes. Here asymptotic flatness will not be 
of much help, because these modes are to be defined by boundary conditions at the past 
horizon. With the lack of spherical symmetry and staticity, one still needs to figure out how 
to make a unique mode decomposition and to obtain a “natural” vacuum state. [^] 

The implementation of this method in the non-spherical case is of course technically 
more challenging (even for a non-BH background), because now the mode functions 
which comprise F (c o, r, 9, ip) are to be obtained by numerically solving PDEs rather than 
just ODEs. 


IV. CALCULATION OF U 2 ) IN SCHWARZSCHILD 

w / ren 


Using the method presented in the last section, we compute (4> 2 ) ren in the exterior region 
of Schwarzschild spacetime, in the Boulware vacuum state, for a minimally coupled massless 
scalar field. The Schwarzschild metric is 


, 9 , 2 M\ , 7 ( 2 M 

ds 2 = — f 1 -) dt 2 + 1- 


-i 


dr z + r 2 dil 2 


(4.1) 


In th is metric the quantum field (j) (x) can be expanded in the form presented in Eqs. (2.2)- 
(2.3), where (r) is conveniently recast as 


Vu (r) = 


r^/Airu) 


~-A,i (r), 


and ipui (r) satisfies the radial equation 

d 2 ^ui (r) 


dr 2 


= - [uj 2 - Vi (r)] -i/w (r) • 


Henceforth, r* will denote the tortoise coordinate given by 


r* = r + 2 M In 


- 1 


and 


is the effective potential. 


V l (r)= 1 - 


2 M 


2 M 


1(1+ 1) 2 M 


(4.2) 


(4.3) 


(4.4) 


17 Furthermore, with an arbitrary choice of r,9,ip coordinates, and with a corresponding arbitrary con¬ 
struction of the set of “up” modes, we have no guarantee that the resultant “vacuum state” would at all 
be a well-defined Hadamard state. In this regard, we should mention the observation that there is no 
Hadamard state that respects the symmetries of Kerr spacetime and regular everywhere. dug 
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The general solution of Eq. 
we denote (r) and ij>™ p (r). 
taken to be 


( 4.3[ ) (for given u, l ) is spanned by two basic solutions, that 
The boundary conditions for these two basic solutions are 


Cl (0 


rue-™*, 

+ Peal 


r* -)• — oo 
r* ->• oo 


V’, 


lip 

cj / 


r = 


+ Peal e" 


r* ->■ -oo 
r* —)• oo 


(4.5) 


where r^, p^z, and p u i = — p^T^i /represent the transmission and reflection amplitudes. 
These two basic solutions are properly normalized and mutually orthogonal. 

The presence of two independent modes for each culm (as opposed to a single such mode 
in e.g. Minkowski) requires a slight modification of the formalism above: We now have two 
sets of annihilation operators, a™ m and a^f lm —as well as their conjugate operators a'"i m , 


a loim- Correspondingly, in Eq. (2.2), in addition to the integral over lu and the sum over l 


and m, we also have to sum over the separate contributions of the “in” and “up” modes. We 
shall consider here the B oulware vacuum , namely the quantum state annihilated by all the 
operators a^ p lm as well as . Revisiting the analysis of the previous sections, one finds that 
everything remains intact, except that all the equations that involve summation over l,m 
should now also include a summation over the “in” and “up” contributions. Correspondingly, 


Eq. (3.3) is now replaced by 


F ^ r > = E (!*” (r) l z + K m 


1=0 


(4.6) 


but otherwise the results of the previous section, and in particular Eqs. (3. T0|3.11|3.12D , are 
unaffected. 

Since the Schwarzschild metric is a vacuum solution, and since we are dealing with a 


massless field, the counter-term (2.6) now reduces to 

G ds (x, x) = 


87rV 


(4.7) 


a is related to s = t' — t through the proper-time r of the short geodesic connecting 1 to i', 
via cr = —t 2 /2. By conducting a second-order expansion of the geodesic equation we obtain 


o(e) = 


1 - 2 M/r _ 2 M 2 (1 - 2 M/r) ei + f 


-£ — 


24r 4 


and correspondingly 


h G D S (x,x ) 47r 2( 1 _2 M/r) 


£~ 2 + 


M 2 


487r 2 r 4 (1 — 2 M/r) 


+ 0(£). 


(4.8) 


(4.9) 


Comparing this to Eq. (3.6) we find that 

1 


a(r) = 


47 t 2 (1 — 2 M/r) 


, d(r) = 


M 2 


487 r 2 r 4 (1 - 2 M/r) ’ 


c(r) = 0. 


(4.10) 
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Therefore in the Schwarzschild case Eqs. (3. T0|3.11|3.12D reduce to 


( 0 2 ( x ))ren = h Fre a ( uj,r)du-hd (r). 


(4.11) 


with 


F reg (w, r) = F (uj, r) + a (r) w. 


(4.12) 


Summarizing the analytical part of the calculation, (</r (x)) ren is given by Eq. (4.11) along 
with Eqs. (4.12), (4.10) and (4.6). Note that basically this expression for Schwarzschild was 
already obtained by Candelas, HU, but here we also complete the calculation by implement¬ 
ing the numerical part as well (and by doing so, we encounter the oscillations problem and 
address it). 


A. Numerical implementation 


We have numerically solved the radial equation for jJj 1 ^ (r) and (r), using the standard 
MATHEMATICA numerical ODE solver, in the domain 0 < uj < 3, at a set of u j values with 
a uniform separation duj = 1/300. Hereafter, we use units in which the BH mass is M — 1 
(in addition to C — G — 1), hence uj as well as r are dimensionless. 

For a given uj, the contribution of the different l modes to F (c o,r) starts to decay expo¬ 
nentially fast beyond a certain l value, typically of order ~ cur (1 — 2 M/r)~ 1 ^ 2 . (This decay 
may be interpreted as tunneling into the potential barrier.) Correspondingly, for each uj 


value, we truncate the sum (4.6) at an l value where the contribution becomes negligible 


(< 10 10 ). Then from F (uj, r ) we construct the regularized function F reg (uj, r ) according to 
Eq. d4T2b. 


Figure la displays F(uj,r) for r = 6, which we choose here as our representative r 
value. Then Fig. lb displays the regularized function F reg (uj,r), obtained from F(uj,r) 
by removing the linear piece a (r) uj. Clearly, the linear divergence has been removed, but 
there remain oscillations that grow as cu 1 / 2 , which we address below. The origin of these 
oscillations (the aforementioned connecting null geodesics), and the determination of their 
frequencies, are discussed in Appendix [Bj 

To calculate {4> 2 (x)) ren we need the generalized integral of F reg (uj,r ) from uj = 0 to 


infinity, see Eq. (4.11). To this end we define the integral function 


H(uj) = 


^reg (^ ; ^*) • 


(4.13) 


The strict integral of F reg would of course correspond to the limit uj —>• oo of H(uj). Figure 
[2] displays H(u), and makes it clear that this limit does not exist, due to the growing 


oscillations (which were inherited directly from F req ). We therefore have to resort to the 


generalized integral instead, as discussed in Sec. IIB The Abel-summation integral (2.10) 


is well defined in this case. However, since we know the precise frequency of oscillations, it 
is much more convenient and more efficient to employ the self-cancellation integral (which 
is fully consistent with the Abel integral). As it turns out, there are multiple oscillation 
frequencies in F reg (uj ) (see Appendix [Bj, and we cancel each of the four dominant ones by a 
forth-order self-cancellation operation. Adapting the notation of Appendix [A] to the present 
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( a ) The numerically calculated 
F (u>, r = 6) in the Schwarzschild case. 


(b) F reg (ui,r = 6), which is the result of 
subtracting the linearly-diverging piece 
a (r) to from F (ui, r). 


Figure 1 



Figure 2: The integral function H(oj,r = 6) (that is, the integral of F reg from zero 

to oj). 


specific context (integration of F reg ), the desired generalized integral is 


"Oo(sC*) 


F reg (oj)doj = lim 


(4.14) 


where = T* [//(w)], H(uS) is given in Eq. (4.13), and T* denotes the multiple self¬ 


cancellation operation as generally defined in Eq. (A6) and detailed in Eq. (Bl). 

Figure [3] displays the function h/*( oS). It is remarkable that after the oscillations have 
been removed the function converges very quickly. This allows a fairly precise determination 


of the limit cu —> oo of this function, which constitutes the generalized integral in Eq. (4.14) 


We then substitute this integral, as well as d(r ) of Eq. (4.10), in Eq. (4.11). In Figs. 4a 
and 4b we present our results for (4> 2 ) ren as a function of r, and compare them to results 
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H„(tj) for r=6 



Figure 3: The function H^(uj,r = 6) (obtained from H after self-cancellation of the 
oscillations). Notice the quick convergence. 


obtained previously by Anderson [13] using a very different method (analytic extension to 
the Euclidean sector). The differences are typically of order a few parts in 10 3 , consistent 
with the estimated numerical errors. 



r 


( a ) Comparing ( 4 >2 ) ren calculated using 
our new method (the red pluses) to the 
previous calculation by Anderson using the 
Euclidean sector and the WKB expansion. 



r 


(b) Zoom-in on the box in Fig. 
the values close to zero. 


4a 


Figure 4 


V. DISCUSSION 

We presented here a new approach for implementing point-splitting regularization numer¬ 
ically, for the computation of (4> 2 ) ren and ( T a p) ren in various asymptotically-flat spacetimes. 
Our main motivation in developing this approach is to allow systematic investigation of self- 
consistent semiclassical evaporation of BHs. This would require the calculation of ( T a p) 
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in a time-dependent BH background. 

So far we developed two variants of our basic approach (both for a quantum scalar field 
in asymptotically-flat background): (i) f-splitting, applicable to stationary spacetimes, and 
(ii) angular splitting, applicable to spherically-symmetric spacetimes. In this first paper we 
focused on the simplest of the two, the f-splitting variant. In presenting this method, we 
restricted our attention to static spherically-symmetric backgrounds, and to (</> 2 ) ren rather 
than (T a0 ) ren , for the sake of simplicity. But we also described the extension of the method 
to more generic, lion-spherical, stationary asymptotically-flat backgrounds. This extension 
suffices for the non-BH case, e.g. for a rotating star. In the case of a stationary BH which is 
not a spherically-symmetric static one, our method still needs a completion: We still have 
to formulate the construction of an appropriate “up” state at the past horizon, which would 
constitute a physically-meaningful (and properly Hadamard) “vacuum” state. 

We then implemented the f-splitting variant to the specific case of Boulware state in 
Schwarzschild spacetime (for a minimally-coupled massless scalar field). The analytical part 
of the regularization procedure coincides in this case with the one developed by Candelas mi. 
However, here we also implemented the numerical part, which involves the numerical solution 
of the radial equation for the various modes, and the summation/integration over the mode 
contributions. Doing so, we found that the regularized function F reg (oj, r ) — which was 
naively expected to be well-behaved at large u ; — actually suffers from growing oscillations, 
which make the ^-integral non-convergent. As it turns out, this phenomenon has little to do 
with the short-distance behavior of the TPF in the coincidence limit x' —» x. Instead, the 
oscillations originate from divergences of the TPF at remote points x' which are connected 
to x by null geodesics. We used the notion of generalized integral — and particularly the 
pragmatic method of self-cancellation integral — in order to handle these oscillations and to 
carry the desired integration over oj. Doing so, we found excellent agreement with previous 
results obtained by a different method, the Euclidean extension [T3j. 

Putting aside for a moment the ultimate goal of analyzing the time-dependent evaporation 
process, we wish to emphasize that even the simplest version of the method, the f-splitting 
variant presented here, makes it possible to do PS regularization in non-spherical stationary 
spacetimes, e.g. that of a strong-field spinning star (although in the generic stationary BH 
case our construction still needs a completion). This was not possible so far, due to the 
difficulties in conducting high-order WKB expansion in such spacetimes. 

However, to achieve our primary goal of self-consistent semiclassical evaporation, we must 
deal with time-dependent backgrounds. To this end we shall need the angular-splitting 
variant, which is slightly more complicated than f-splitting. We shall describe this variant 
in a separate paper. We already applied the angular-splitting method to (< f> 2 ) ren and also to 
(T a0 ) ren in Boulware state in Schwarzschild, and again we found very good agreement with 
previous calculations [T3J. These results will be presented elsewhere [IB] . 


ACKNOWLEDGMENT 

We are grateful to Paul Anderson and Bob Wald for many interesting and fruitful discus¬ 
sions. In addition AO thanks both Paul Anderson and Bob Wald for kind hospitality in their 
institutions, where many of these discussions took place. In particular the understanding 
of the nature of the oscillations, and their relation to connecting null geodesics, emerged in 
a discussion with Bob Wald during a visit at the University of Chicago. We also wish to 
thank Paul Anderson for sharing his unpublished data with us. 


19 


Appendix A: Multiple self-cancellation 


For a given integrand function h(oj), the self-cancellation integral was defined in Sec. 2.2 


as 


where 


T\ is defined by 


r*oo(sc) 


Tx [/Ml = 


h(uj)duj = lim T\ [H(u >)\, 

id —^OO 
pdd 

H(u>) = / h(x)dx , 

Jo 

f( u, ) + f(u> + A/2) 


(Al) 


and A is the oscillation’s wavelength. 

Consider next the situation in which h(uj ) contains oscillations in two different frequencies 
(say, with u -independent amplitudes). Then H{uj ) will take the asymptotic form 


H(u) ~ Aie lulLl + A 2 e lU}L2 + const 


(A2) 


with two different wave-lengths Ai i2 = 27 t/L 1>2 . If we apply a self-cancellation with respect 
to (say) Ai, then we will be left with a function T\ 1 [H(uj)} in which the const is of course 
unaffected, and the term oc e luLl has entirely been annihilated. The other oscillatory term 
oc e lbjL2 is still present in T Al [//(u;)], although its amplitude will decrease by the factor 
cos [(L 2 /Li)tt/2], as one can easily verify. The self-cancellation integral (associated to TxJ 
is thus still iron-convergent. A second self-cancellation, this time with respect to X 2 , will 
yield a convergent self-cancellation integral, which may be formulated as 


lim T A2 T Al [H(u)} 


The order of the two operations Tx, and T\ 2 is unimportant, as one can easily verify. With 
both orderings, the double self-cancellation integral yields the “const” in Eq. (A2). [^] 

This process is straightforwardly generalized to functions H(uj) with any number k of 
oscillation frequencies that need to be annihilated. To handle this situation, we define the 
multiple self-cancellation operator 


T* = T\ 1 T\ 2 ...T\ k . 

Then the multiple self-cancellation integral may be expressed as 

/»oo(sc*) 


where 


h(u)du = lim 

CJ^-OO 

H*(cj) = T* [H(u>)\ . 


(A3) 

(A4) 

(A5) 


18 To avoid confusion, we point out that an application of self-cancellation operation T A /, but with a “mis¬ 
taken” wavelength parameter A' that differs from the true oscillation wavelength A , does not “spoil” the 
generalized integral in any way: The mismatch in A does not lead to any new oscillatory terms, nor does it 
modify the value of the generalized integral. The effect of the mismatch in A is merely to limit the efficiency 
of the self-cancellation operation: It decreases the oscillation’s amplitude by the factor cos [(A'/A) 7 t/ 2] 
instead of fully annihilating it. 
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As was demonstrated in Sec. IIB if H{u) contains an oscillation with fixed amplitude, 


then the T\ operation fully nullifies this oscillation. (This was demonstrated there for the case 
h(oo) = sin(o;L), but it equally applies to the more general cases h(oo) = e ±luL .) However, if 
the oscillation’s amplitude varies with c o, the cancellation is not complete. Consider the case 
H{ui) = g(u)e lujL where g(oo) is some slowly-varying function — namely, a function whose 
typical length of variation £ becomes A at large oo. To be more specihc, let us further 
assume that £ diverges as oo — y oo. [j^] Then one can easily show that 

at large oo. In particular, if H(u) = oo p e lujL then 

T x [H(u)} oc pu p ~ 1 e iojL . 


Consider now the case H(oo) ~ oo p e luL for some 1 < p < 2. Then T\ [H(oo)\ still does not 
converge as u —> oo. Nevertheless, T\ (T\ [H(oo)}) is oc oo p ~ 2 e llJjL and hence it converges. This 
illustrates that in certain circumstances one may want to apply the same self-cancellation 
operation several times, say n times, an operation which we shall denote as (T\) 

oo p e luL one finds that 


Quite generally, in the case H(uo) 
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(T x ) n [H(u})\ oc u- 


,P~n e iuL 


Therefore, for this class of H(uo) functions, the convergence criterion for the (single¬ 
frequency) multiple self-cancellation integral is simple: It converges if and only if n > p. 

In the case 0 < p < 1 a single self-cancellation operation would be sufficient for achieving 
convergence. However, one would then be left with a slowly-decaying amplitude oc oo p-1 . To 
speed the large-w convergence, one may repeat T\ several times. This is in fact the situation 
in our specihc problem, where h(u) oc (n 1,/2 e* aji and hence H{uo) oc u l ^ 2 e lu:L + const: The 
single self-cancellation integral converges, but rather slowly, with amplitude oc a; -1 / 2 , and we 
therefor repeat T\ several times to achieve faster convergence (which is crucial for numerical 
implementation). 

If there are several different oscillation frequencies that need to be annihilated, one can 
freely choose how many times to repeat the T x operation for each frequency. We therefore 


generalize the above expression (A3) for the multiple self-cancellation operator: 


T* = (T Xl r (T A2 r...(T A j 


\nk 


(A6) 


The multiple self-cancellation integral (sc*) is thus obtained by using this T* operator in 
Eqs. ( A4|X5 ). 

Finally we point out that this multiple self-cancellation integral is fully consistent with 
the Abel-summation integral — in the same sense discussed in Sec, IIB (concerning the single 
self-cancellation integral). 


19 This includes for example all powers oj p , because in this class £ oc oo. (But it also includes much more 
general classes of functions, e.g. w p (lnu;) 9 for any p, q.) 

20 The exception is the case of natural p with n > p, in which ( T\) n [-ff(w)] strictly vanishes. 
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Appendix B: The oscillations in F reg (u,r ) 


As can be clearly seen in Fig. lb, the function F reg (a;,r) admits growing oscillations in 
oj. These oscillations are of course inherited directly from F (c u, r), as can be seen in Fig. la 
(although in the latter the oscillations are overshadowed by the linearly-growing term). In 
this Appendix we shall discuss the origin of these oscillations, their frequencies, and their 
amplitudes. Then at the end we shall specify the self-cancellation operator that we apply 
in order to practically remove these oscillations. 


1. Origin and Nature of oscillations 


Owing to spherical symmetry and staticity, in t-splitting the TPF may depend only on r 
and on e — t' — t. Let us then introduce the abbreviated notation for the TPF 

P(e, r) = (0 (t, r, d, <p)(p(t + e, r, 0, <p)). 

Equation (3.4) actually tells us that F (ca, r) is the Fourier transform of P(e, r). If P(e, r ) 
were regular and smooth for all e, then its Fourier transform F (ca, r) would decay quickly 
at large u, faster than any power of the latter. The undamped oscillations in F (ca, r) must 
therefore indicate some irregularity in P(e,r). We still need to understand the nature of 
this singularity, and its location on the £ axis. 

The singularity of the TPF at £ —> 0 indeed leads to a divergence of F(uj,r) at large 
oj (this is the linearly-growing term shown in Fig. la) — but not to oscillations; And this 
linear singularity is no longer present in F reg ( oj,r ). The oscillations in Fig. lb must then 


indicate another singularity, located at some point e = e s 0. To illustrate this, consider for 
example the delta-function case: The Fourier transform of 8{e — e s ) is e~ lU£a . It is oscillatory 
if and only if e s 0. The frequency of the cn-oscillations is just e s — namely, it directly tells 
us the distance of the singularity from the point e = 0. 

One more example is the singularity \e — £ s \~^, whose transform is oc co^ 1 e~ lU}£a . But this 
phenomenon is of course more general: If a certain function P(e) admits a Fourier transform 
F (ca), then the transform of P{e — e s ) is F (cj) e~ lU£s . 

The oscillations seen in e.g. Fig. lb have a certain “cn-wavelength” Ai ~ 0.17 (in units 


in which the BH mass is M — 1). They must therefore correspond to a singularity in the 
TPF, located at a distance e s = t! — t = 2n/\i ps 37 from point x. What is the nature 
of this non-local singularity of the TPF? As already pointed out in Sec. |IIB[ the function 


(0 (x) 0 (x')) admits a singularity whenever a null geodesic exists which connects x and x'. 
In the present context of Schwarzschild background and ^splitting, we are dealing here with 
a null geodesic which emanates from a certain spatial point (r,9,ip), makes a round trip 
around the BH, and then returns to that same spatial point, but obviously with a certain 
delay in t — which should correspond to the shift parameter e s . 


2. Spectrum of oscillations 

In fact there is not only one but an infinite, discrete set of such connecting null geodesics 
(for each r). This is because a null geodesic emanating from a spatial point (■ r,9,(f) can 

21 In this context we should regard F (ui,r) as a function that vanishes for all oj < 0. Note also that in 
the present context one should not think of e as a small parameter: Instead, it is allowed to take all real 
values. 
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make any integer number of revolutions around the BH before returning to that spatial 
point. Therefore, the TPF will actually admit an infinite number of singular points at a 
discrete set of values t' — t — e n , one for each integer n. 

Correspondingly, there will be a discrete spectrum of oscillation modes in F reg (oj, r), with 
(r-dependent) oj -frequencies e n and corresponding wavelengths X n = 2 tt/ e n . We point out, 
however, that the dominant mode is always n = 1, and the oscillation’s amplitude quickly 
decays with n (see next subsection). 

To perform the self-cancellation of oscillations, we shall need to know the spectrum of 
frequencies £ n , at any desired r value. This requires integration of the null geodesic equation, 
to find the connecting null geodesics. For r = 3 M the situation is especially simple, because 
in that case the connecting null geodesic is circular. One then finds that e n = (2ny/27) n. 
At other r values the spectrum becomes more complicated, and needs to be calculated nu¬ 
merically. The null orbits in Schwarzschild are characterized by a single constant of motion, 
namely the angular momentum per unit energy (or the “impact parameter”). The calcu¬ 
lation of e n involves (i) numerical integration of the null geodesic equation (for prescribed 
values of that constant of motion), and (ii) using the Newton-Raphson method for adjusting 
this constant, in order to find the connecting geodesics, which return to the original spatial 
point after n rounds. The required parameters £ n are the f-duration of these connecting null 
geodesics. Overall, this is an easy numerical procedure. 

In our representative case r = 6 M, the first few oscillation frequencies are found to be 


ei ~ 37.50 , e 2 - 70.17, e 3 ^ 102.8 , £4 ~ 135.5 


In Fig. lb we predominantly see the basic oscillation n — 1, and the frequency agrees very 
well with this value of £4. We can also notice the residual effect of the n — 2 oscillation, 
which causes the small distortion (i.e. small deviation from the smooth oc cu 1//2 envelope) 
in the pattern of peaks of F reg (c o,r), seen in Fig. lb After the basic oscillation n — 1 is 
removed by self-cancellation, the next one (n = 2) dominates and becomes very clear, and 
again, its frequency is found to agree very well with the above value of £ 2 . By this procedure 
it is possible to expose a few more modes, and to confirm their agreement with the above 
£ n values, that were obtained from the connecting null geodesics. 

As was already mentioned above, for r = 3 M the frequencies e n form an arithmetic 
sequence with a common difference Aeo = 27 ta/ 27 ~ 32.65. For other r values this is no 
longer the case. Still, the difference £ n+ \ — £ n quickly approach the standard spacing Ae 0 - 
For example, from the four £ n values specified above for the r = 6 M case, one sees that 
£2 — £1 ~ 32.67, and then e 3 — e 2 ~ £4 — £3 ~ 32.65. This is simply because for any r, a 
connecting null geodesic with large n makes most of the revolutions around the BH along 
an orbit very close to the circle r = 3 M. 


3. Amplitude of oscillations 

The numerical data indicate that the oscillation’s amplitude grows as a; 1 / 2 . This in turn 
implies that the divergence of the TPF at £ —> £ s should be oc |e — e s | -3 / 2 . We shall not 
address this issue here in detail, but we point out that this —3/2 power is just what one would 
expect from simple arguments. To this end one has to recall that since the background is 
spherically symmetric, and since the splitting is only in the t direction (implying that O' = 0 
and ip' = <£>), whenever x' is located on a null geodesic emanating from x, it is placed exactly 
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at a caustic point of that null geodesic. Qualitative arguments suggest that on crossing such 
a caustic point, the TPF should indeed diverge as \e — e s Q 3//2 . But the discussion of this 
issue is far beyond our present scope. 

The n > 1 oscillations, too, grow as a; 1 ^ 2 at large u (the numerics confirms this, at least 
for the first few n values). This is for the same reason as that described above for n — 1. 
We thus express the large-cu amplitudes of the various e n modes as ~ A n cu 1 / 2 , where A n is 
a set of (r-dependent) amplitude parameters. 

It is important to explore how A n behaves with increasing n, in order to control the 
possible effect of the infinite number of oscillating modes. The simplest case to analyze 
is again r = 3 M, because the connecting null geodesic is then the circular geodesic at 
r = 3 M. Simple analytical considerations suggest that in this case A n should form an 
almost-exact geometric sequence with A n /A n+ i = e™ ~ 23.14. [^] Our numerical results for 
F reg (c u, r = 3 M) allows reliable evaluation of the first three amplitudes, and the calculated 
ratios A1/A2 and A^j'A3 agree very well with e 7r , to about one part in 10 3 . (For the n > 3 
modes the oscillations are too weak to reliably measure their A n .) 

For other r values the situation is more complicated, and we do not expect to find such 
a well-approximated geometric sequence; Yet, we still expect that as n increases, A n /A n+ 1 
should quickly approach the above “canonical” value e n . The reason is that, for large n, 
the connecting null geodesic makes most of the n revolutions around the BH along an orbit 
very close to the circle r = 3 M. Hence, the decrease in the Van-Vleck determinant at each 
revolution is approximately the same as in the analogous r = 3 M case (an approximation 
that ever improves with increasing n). At very small n, however, the ratio between two 
successive amplitudes may slightly differ from eA. We numerically find that the ratio A 1 /A 2 
ranges from e A ~ 23.1 at r = 3 M to ~ 23.7 at r = 9 M. 

Overall, at least in the range 3 M < r < 9 M that we have numerically explored, the 
numerical data as well as the theoretical considerations are all consistent with an almost- 
geometric sequence (even for small n), with A n /A n+ 1 ranging between 23 and 24. In turn 
this implies that for practical computation of the generalized integral of F reg (co, r), we shall 
have to cancel the first few n modes, but the contribution of large-n modes may be neglected. 


4. Self-cancellation operator 


We self-cancel the first four frequencies n = 1...4. The higher modes n > 4 are too weak 
to notice. ^ The values of the frequencies e n are numerically obtained from the connecting 
null geodesics, for any desired r, as explained above. The corresponding wavelengths are 
then given by X n = ‘I'nje n . For each frequency, we apply a fourth-order self-cancellation. 


Thus, in the terminology of Eq. 

T, 


(A6), our actual multiple self-cancellation operator is 
^(T Al ) 4 (T A2 ) 4 (T A3 ) 4 (T A4 ) 4 . 


(Bl) 


22 A simple (though still unproved) analytical argument, based on evaluating the Van-Vleck determi¬ 
nant along the r = 3 M geodesic, suggests that in this case A n should be exactly proportional to 
[2sinh(27rn)]^ 1,/2 . This expression deviates from the geometric sequence e _7ra by a tiny relative amount 
= e~ im /2. This deviation is smaller than one part in 10 5 even for n = 1, too small to be detected by our 
numerics, but nevertheless our numerical results are fully consistent with that expression. 

23 We point out that although the amplitude formally diverges at large cu for any n, because after all we 
only integrate along a finite interval of to [which is in turn allowed due to the fast convergence of f/*(w)[, 
and because of the fast decay of A n with n, the large-n oscillation terms do not have any noticeable effect 
on the integral. 
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